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This report sunmar izes the results of the work performed by Advanced 
Computer Concepts. Inc. under Contract NAS5-31218 with the National 
Aeronautics and Space Administration Goddard Space Flight Center 
titled ”GHRS Science Verification Data Analysis.” 


1. CALIBRATION SUPPORT FOR TUT GHRS. 

The majority of the data analysis we performed to support the Orbital 

Ve r i f 1 c a t ; on ( (>Y i and Science Ye r i f i c a t i or. ' SY 1 o f t he Godda r d High 

Resolution Spectrograph (GHRS' was in the areas ol the Digicon 

detector‘s r e r f o rm^ nee and stability. wavelength calibration, and 

£ e oma £ n e t i c induced i ma g e mo t i o n . The fell ov : n g sections briefly 
describe the results of the analyses. Detailed results can be found 
in the a t : a w hme n t s . 


1 . ] D i g i c c n P e r f o rma nee and Stability. 


The geometric transformation which relates the Dig icon's X and V 
deflection* and the diode position to a position on the photocathode 
is modeled by the equations ir. Attachment 1. We computed the 
c oe 1 f i c i e n t s of this model using both pre-launch calibration data 
taken in 19S4 and post-launch data taken during 0\ and S\. Our 
analysis showed no significant changes with time (Attachment 1 and 2). 
This indicates that any aging effects in the Dig icon’s permanent 
ma g n e t is m i n i ma 1 . 


A detector stability mi n i f unc t i ona 1 wa s run 10 times for both GHRS 
detectors during the period of November 12, 1990 through June 9, 1991. 
Our results show no changes in the focus or the ge ome trie properties 
of the detectors over the 7 month period. The deflection step size 
remained stable to within 1 percent and the dark count remained at 
approximately 0.01 counts second diode for both detectors. The only 
significant change was a decrease in the observed flat field lamp 
count rate. The count rate is decreasing at a rate of approximately 7 
percent per year. It has not yet been determined whether this 
decrease is a result of a detector sensitivity loss or a decrease in 
the light output of the flat field lamps. Simular decreases were 
found in the data from the wavelength calibration lamps which would 
indicate that a loss of detector sensitivity is the most likely 
explanation. More detailed discussion of the detector stability can 
be found in Attachment 2. 


1.2 Geomagnetic Induced Image Motion. 

We have found that the GHRS detectors are susceptible to 
geomagne t i c a 1 1 v induced image motion from ana>s:s of OX SY spectral 
calibration lamp data and target acquisition deflection calibration 
data. We have measured the motion due to the Forth’s magnetic tield 
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and found that the mo tion at the diode array is a p p r o x 1 ma t e 1 y 10 
microns (0.2 d i ode s -'Gaus s ) for detector 1 and 50 microns (1.0 
diode Gauss) for detector 2. This motion can cause errors in target 
acquisition. errors in measured wavelengths, and can result in a loss 
of spectral resolution. The impact of the magnetically induced motion 
can be minimized by limiting spectral observations to less than 5 
minutes and by using the double located mode of target acquisition. 
Detailed results of the magnetic motion problem can be found in 
a 1 1 a c lime n t 3 . 

1.3 ( )Y SY Wa \ e 1 e n g t h Calibration. 

We have analyzed the GHKS OY SY spectra! calibrator 1 ..mn observations 
and have found that an average dispersion coefficient table with 
1 : n e a r the rma 1 mo t i o n mo dels can be used to c ompu te wavelengths for 
science observations with the target in the small science aperture 
(SSA) to better than one photocathode sample unit '.50 microns). We 

have found that the ma \ o r i t y of the spectral image mo t i o n can be 

modeled as a linear function of the temperature. Selection of the 
best Thermistor to use in the model varies with grating mode. 

Gratings 3 and 4, with a motion range of 300 microns, are most 

susceptible to the rma 1 mo tion. J n a d d i t i o n to the the rma 1 mo t ion, we 
have observed plate scale variations on the order ol 0.1 percent which 
ma y a Iso be the rma 1 1 y induced. De tails of this analysis can be found 
in a t t a c hme n t 4 . 

We have ma d e me a s u r erne nts of the spectral offsets be twe en the large 
iLSA) and small science apertures for GHRS gratings 2, 3, and 4. In 
all cases the measured offsets were less than one diode and the 
average for each grating was less than 1 ' 2 d i d o e . There was 
insufficient data to quantify changes in the offsets for varing 
carrousel and photocathode sample positions. However, our results for 
grating 2 indicate that the offset does vary with carrousel position. 
More detailed analysis of the LSA to SSA offsets can be found in 
a r t a c hmen t 5 . 


2 . 0 SOFTWARE DEXTLOPMENT 

Another major task in our support of GHKS OY SY data analysis was to 
provide software developemnt and maintai nance support for the GHKS 
data management and analysis system. This system is installed on the 
GHRS Yax at NASA/Goddard Space Flight Center running under the VMS 
operating system. A significant portion of the s o * t wa r e vva s also 
converted to run under both UNIX and MS-DOS operating systems. The 
s o f twa re i s wr i t ten a 1 mo s t exclusively in Interactive Da t a Language 
( IDL) . 

We have continued the conversion of the GHKS sol twa re : r oni IDL version 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Page 3 


1 to 1DL version 2. The major reason for this conversion is that IDL 
version 1 only runs on a VAX WiS machine. Conversion to version 2 
allows use of the software under \MS . CM\ and .\1$-IX)S machines. 
Although only minor changes in most of the software was required to 
convert to IDL version 2 a significant amount of effort was required 
to convert internal binary data sets fr om one machine to another. We 
have written routines to convert the three mo st popular data formats 
used by the GHRS software system from the VAX to other machines. 
These formats include the Sc ience Data Analysis System t SDAS ) image 
and table formats and the GHRS data base format. 


o . < GI K S CAL 1 BRAT ION RL FERENCF F i LH S . 

W'c 1 v rma t t e d the f o • i owi ng calibration files c ompu ted b> the Gl IK S t e am 
according to the specifications in the Space Telescope Science 
Institute (ST Sc I > document ’'Hubble Space Telescope Spectrograph 
Ca 1 i b r a t ; o i: ITo t a Specifications.” We delivered the t c the ST S c 1 for 

use the the routine data processing of GlfRS data. 

FILL Description Delivery Date 

DC__002.TAB Average OY SY wavelength dispersion June 10. 1991 

coefficients. They do not include 
a tiie rma] motion model. 

LA_00 1 .TAB Incidence angle corrections for June 10, 1991 

c omp u t i n g the wa velength offset 
between the small science aperture 
and the spectral calibration lamp 
apertures . 

S\IAP_002 . TAB Average photocathode sample mapping June 10, 1991 

function coefficients. 

YG_G1_01 .HH* Files containing the vignetting June 10, 1991 

VG__G2_02 .HH* curves for the first order gratings. 

VG_G3_01 . HH* 

VG_G4_01 -HH* 

VG_G 5__0 1 .HID 

DIO lD1_1 .HH* 

DIO_P2_l .HH* 

SF_G1_01 .HH* 

SF_G2_01 .HH* 

SF_G3_01 - HH* 

SF_G4__01 - HH* 

SF_G5_01 -HH* 
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Contains the diode response computed June 10. 1991 
from internal flat field lamp 
observat ions . 

Contains the sensitivity curves for June 10. 1991 

each of' the first order gratings. 



June 10. 1991 


WT_G1_01 .HH* 

\VF_G2 01 .HH* 

\VT_G3_01 .HH* 

WF_G4 01 .HH* 

\YF_G5_01 . HH* 

1 X' ('04 . TAB 


. TAB 


i.A_""l .'JAB 


S\1AP_003 . TAB 


YG_G1_02 .HH* 

\G_G2 03 .HH* 

\G_G3_02 .HH’ 

YG_G4 02 .HH’ 

\G_G5_02 .HH* 

\V1_G1_01 .HH* 
\VF_G2_01 .HH* 
\\T_G 3_0 1 .HH’ 
\M_G4_01 .HH* 
\\F_G5_0l .HJJ* 
SF_G1_01 .HH* 
SF_G2_01 .HH* 
SF_G3_01 .HH* 
SF_G4_01 .HH* 
SF_G5_01 .HH* 

\VF_F1A_01 .HH* 
SF_XA_01 .HH* 
\VF_EB__01 .HH* 
\\F __EB_0 1 .HH* 

CXTTA_2 . TAB 
CTJR 9_2 . T.AB 


] A_0(>2 . T.AB 


Contains the wavelength vectors for 
the above sensitivity curves. 


N ew 
coe 
f o r 

average wavelength dispersion 
fficients which include provisions 
a thermal motion model. 

Sept . 

6 . 

1991 

('ce 
mo- 1 

fficients for the rma 3 spectre;] 
ion. 

Sept . 

6 . 

1 9 9 1 

Incidence a n g 1 e coefficients Tabulated 
with new grating names. Actual 
coefficients are unchanged. 

Sept . 

o . 

1 991 

1 n 
c o e 

ai SY average sample mapping 
fficients. 

Sept . 

6 . 

1 991 

L pd 

ate of precious files. These 

Sept . 

6 . 

1991 


files are constructed so that 
interpolation between curves is 
done using the photocathode line 
position instead of carrousel position 

Redeli eery of the first order Sept. 6. *991 

grating sensitivity curves using 
the n ew g r a t i n g n ame s . 


Sensitivity curves for the echelle Sept. 6. 1991 

g r a t i ng mod e s . 


Tables containing improved echelle B Oct. 30. 1991 

ripple coefficients computed from 
prelaunch miniarc 1 amp data. 

Contains the incidence angle Oct. 30. 1991 

coefficients for gratings 2. 3 and 4 
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large science aperture data. 

CXT\3_1 .TAB Table of detector parameters which Dec. 1. 1991 

contain the order of the new 
background fitting polynomial. 

C.XKB . TAB Be he lie scattered light coefficient Dec. 1, 1991 

table . 

4 .m restoration or iras large Magellanic clo vn data. 

We h a \ e continued the res: o ration of t hi e I RA S Add ; t i : ' r. u 1 ( )b servat : o n v 

( AO ' t 1 : the Large Ma gel 1 a r. i c Cloud 1 L\1C ) using t h e ^ * o c k e d iterative 

resii'i a : ; r. a 1 r . r * t hm described in a t t a c hme n • c . Our restoration 
a ] g c- r j t hm wv r k s with the raw data scans \\ i t h r u : r eg r i dd i ne or 

i n t e r p i 1 a t ; n i : ) e data on an e q u a ! 1 v spaced image grid. To do this, 
we must us r a different point spread function for each scan angle, 
direction c. n d d e : e c t o r . 

Results oi the restorations of the 30 Do r a d u s region o : the LML are 

shown in figure 1 for each cf the 4 IRAS wavelength ba nds . The images 
on tiie 1 e f t s h ow the results of regridding the raw data without 
attempting any deconvolution. The images on the right show the 

deconvolution results. figure 2 sh ows the deconvolution of the entire 
data set for the 60 micron wavelength band. 


r 


" £ JS 

quality 



FIGURE 2 

Restoration of the 60 micron IRAS AO data 
of the Large Magellanic Clouds 
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FIGURE 1 

Restoration of the 30 Doradus 
Region of the IMC 


12 microns 


25 microns 


60 microns 


100 microns 


Original Restored 
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ATTACHMENT I 


GHRS PHOTOCATHODE MAPPING FUNCTION 


D.J. Lindler 
11 July 1991 


ABSTRACT 

The post-launch GHRS photocathode mapping function which relates 
digicon deflections and diode positions to photocathode line and 
sample positions show no significant change from pre-launch 
computations. This indicates that any aging effects in the permanent 
magnet or the digicon is minimal. Observations over an 8 month period 
after launch also show no significant changes with time. This report 
describes the post-launch calibration. Tables of the resulting 
coefficients are presented. 


1. I NTRCOXT ION 

The GHRS phocathode line and sample mapping functions which relate 
digicon deflection and diode number to photocathode line and sample 
positions are given by: 

L «= LO + A*dy 
S ~ SO + B*dx + E*D 

whe re , 

L is the line postion in 50 micron sample units 

defined so that the upper photocathode mask edge 
is at pos i t ion 0.0. 

S is the sample position in 50 micron sample units 
defined so that the left photocathode mask edge 
i s at pos i t ion 0.0. 

dy is the y-deflection minus 2048. 

dx is the x-def lect ion minus 2048. 

D is the diode position starting at 0 for the first 
diode on the main diode array. 

LO and A are cofficients which vary by detector. 

SO, B, and E are coefficients which vary by detector 
and y-deflection. 

LO and SO give the photocathode line and sample position which the 
first diode observes at null x and y defelctions (2048, 2048). A and 
B give the deflection step size as a fraction of 50 micron sample 
units. E gives the separation of the position on the photocathode 
viewed by neighboring diodes. 
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These coefficients are computed using observations of the four 
photocathode mask edges illuminated by a flat field lamp. LO and A 
are computed from the y-def lec t ions which center the diode array on 
the top (YTOP) and bottom (YBOTTCM) photocathode mask edge. 

A - 430.0 / (YBOTTCM - YTOP) 

LO = (2048 - YTOP) * A 

The coefficients SO* B, and E are computed from observations of the 
left and right photocathode masks edges observed at various 
x-def lect ions . The diode positions of the edge are tabulated versus 
x-def lect ion. The sample mapping coefficients are then computed from 
the table values by a least squares fit. 


2. DATA REDUCTION 

The line mapping function coefficients were computed using the routine 
MAPLCAL on the observations in table 1. MAPLCAL was called using all 
default parameters by; 

LMAPCAL, [ID1, ID2], 0* TABLE 

where ID1 and ID2 are the observation entry numbers for the top and 
bottom edge observations. The coefficients are written into an SDAS 
table file specfied by TABLE. LMAPCAL was run on each pair of top and 
bottom edge locations listed in table 1. For example, the first pair 
of observations were calibrated using; 

LMAPCAL, [34, 35], 0, ’tablet’ 

The sample mapping function coefficients were computed using the 
routine MAPSCAL on the observations in table 2. MAPSCAL was called 
using all default parameters by; 

SMAPCAL, ID, 0, TABLE 

where ID is the observation entry number and TABLE is the output SDAS 
table file name where the output coefficients are tabulated versus 
y-def lec t ion . 


3. RESULTS 

Table 3 shows the tabulated line mapping function for each pair of top 
and bottom edge scans and table 4 shows an average of the results. 
Small changes in LO can be attributed to statistical errors in the 
measurements and to geomagnetic image motion. The deflection step 
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size vaule, A, is stable to 1 percent. Small changes in this value 
can be attributed to statistical errors in the measurement s , small 
fluctuations in the photocathode high voltage, differential magnetic 
field vectors between the top and bottom edge observations, and 
possible small thermal effects. Table 4 also shows the average value 
of the coefficients computed during the GHRS phase 6 calibration in 
1984. The small changes since 1984 may be the result of adjustments 
to the photocathode high voltage to improve detector focus and digicon 
aging over seven years. 

Table 5 shows the sample mapping function coefficients for each of the 
observations in table 2 tabulated versus detector and y-def 1 ec t i on . 
Table 6 shows the averages over the observations. The deflection step 
size, B, is within 1 percent of the nominal 1/8 diode value (0.125) 
needed to insure that combaddition does not degrade spectral 
resolution. Table 6 also shows that the deflection step size is 
stable to within approximately one percent. Changes in the 
coefficients can be attributed to statistical errors in the 
measurements, small fluctuations in the photocathode high voltage, 
geomagnetic image motion, and possible small thermal effects. Again 
the small deviations from the 1984 coefficients shown in table 7 can 
be attributed to adjustments in the photocathode voltage and possible 
d i g i con aging. 



TABLE 1 

Line Mapping Function Observations 


ENTRY ROOTNAME 

DET 

START_TIME 

Photocathode 

Edge 

34 

Z06H0609T 

1 

16- JON- 1990 

06:52 

Top 

35 

Z06H060AT 

1 

16-JUN-1990 

07:02 

Bottom 

135 

Z06H5609T 

1 

24- JUL- 1990 

10:44 

Top 

136 

Z06H560AT 

1 

24-JUL-1990 

10:54 

Bottom 

565 

Z06H6609T 

1 

10-SEP- 1990 

12:01 

Top 

566 

Z06H660AT 

1 

10- SEP- 1 990 

12:11 

Bottom 

907 

Z0E90N03T 

1 

21 -NOV- 1990 

01:32 

Top 

908 

Z0E90N04T 

1 

21 -NOV-1990 

01 : 41 

Bottom 

61 

Z06H0709T 

2 

18-JUN-1990 

06:27 

Top 

62 

Z06H070AT 

2 

18- JUN-1990 

06:37 

Bottom 

239 

Z06H5709T 

2 

24 - JUL- 1 990 

14:09 

Top 

240 

Z06H570AT 

2 

24- JUL-1990 

14:18 

Bottom 

575 

Z06H6709T 

2 

10-SEP- 1990 

15:26 

Top 

576 

Z06H670AT 

2 

10-SEP- 1 990 

15:35 

Bottom 


TABLE 2 

Sample Mapping Function Observations 


ENTRY ROOTNAME DETECTOR START_TIME 


33 

Z06H0608T 

1 

16- JUN-1990 

06:44 

276 

Z06H5608T 

1 

24 - JUL-1990 

10:36 

564 

Z06H6608T 

1 

10-SEP- 1990 

11:53 

906 

Z0E90N02T 

1 

21 -NOV- 1 990 

01:24 

60 

Z06H0708T 

2 

18- JUN-1990 

06:19 

238 

Z06H5708T 

2 

24- JUL-1990 

14:00 

574 

Z06H6708T 

2 

10-SEP- 1990 

15:17 



TABLE 3 

GHRS Line Mapping Function Coefficients 


De t e c to 

r date 

ID1 

LO 

A 

YTOP 

YBOTTCM 

i 

16- JUN-1990 

34 

214.17 

0.1262 

350.5 

3758.6 

i 

24-JUL-1990 

135 

214.27 

0.1262 

349.7 

3757.9 

i 

10-SEP- 1 990 

565 

214.14 

0.1261 

349.9 

3759.8 

i 

21 -NOV- 1990 

907 

214.24 

0.1263 

351.3 

3756.7 

2 

18-JUN-1990 

61 

225.13 

0.1240 

232.7 

3700.0 

2 

24-JUL-1990 

239 

225.10 

0.1239 

231.2 

3701.8 

2 

10-SEP- 1990 

575 

225.31 

0.1241 

232.0 

3697.8 


TABLE 4 



Average Line Mappi 

ng Funct 

i on Coe f f i c i en t s 


Detector 

LO 

RMS(LO) 

A 

RMS(A) 

L0( 1984) 

A( 1984) 

1 

2 

214.21 

225.18 

0.06 

0.12 

0.1262 

0.1240 

0.0001 

0.0001 

213.56 

224.84 

0.1267 

0.1245 


TABLE 5 

GHRS Sample Mapping Function Coefficients 


Detector date 

ID 

YDEF 

SO 

B 

E 

i 

16- JUN- 1990 

33 

500 

30.17 

0.1260 

1 .00161 

i 

24-JUL-1990 

276 

500 

29.99 

0.1259 

1 .00179 

i 

10-SEP-1990 

564 

500 

30.16 

0.1259 

1 .00166 

i 

21 -NOV- 1990 

906 

500 

30.20 

0.1260 

1 .00168 

i 

16- JUN- 1990 

33 

1000 

30.68 

0.1259 

1 .00192 

i 

24 - JUL- 1 990 

276 

1000 

30.52 

0.1257 

1 .00209 

i 

10-SEP- 1990 

564 

1000 

30.69 

0.1258 

1 .00202 

i 

21 -NOV- 1990 

906 

1000 

30.71 

0.1259 

1 .00211 

i 

16- JUN- 1990 

33 

1500 

31.04 

0.1256 

1 .00218 

i 

24 - JUL- 1 990 

276 

1500 

31.09 

0.1257 

1 .00182 

i 

10-SEP- 1990 

564 

1500 

31 .20 

0.1258 

1 .00187 

i 

21 -NOV- 1990 

906 

1500 

31 . 15 

0.1258 

1 .00211 

i 

16- JUN- 1990 

33 

2000 

31.41 

0.1257 

1 .00224 

i 

24- JUL- 1990 

276 

2000 

31.39 

0.1255 

1 .00231 

i 

10- SEP- 1 990 

564 

2000 

31.51 

0.1256 

1 .00227 

i 

21 -NOV- 1990 

906 

2000 

31.50 

0.1257 

1 .00229 

i 

16- JUN-1990 

33 

2500 

32.05 

0.1258 

1 .00169 

i 

24 - JUL- 1 990 

276 

2500 

31.99 

0.1256 

1 .00195 

i 

10-SEP- 1990 

564 

2500 

32.10 

0.1257 

1 .00188 

i 

21 -NOV- 1990 

906 

2500 

32.04 

0.1257 

1 .00204 

i 

16- JUN- 1990 

33 

3000 

32.53 

0.1255 

1 .00228 

i 

24- JUL- 1 990 

276 

3000 

32.56 

0.1254 

1 .00236 

i 

10-SEP- 1990 

564 

3000 

32.67 

0.1255 

1 .00226 

i 

21 -NOV- 1990 

906 

3000 

32.57 

0.1256 

1 .00245 

2 

18- JUN- 1990 

60 

1000 

28.95 

0.1249 

1 .00445 

2 

24- JUL- 1990 

238 

1000 

29.02 

0.1249 

1 .00443 

2 

10-SEP- 1990 

574 

1000 

28.76 

0.1248 

1 .00457 

2 

18- JUN- 1990 

60 

1500 

28.98 

0.1246 

1 .00295 

2 

24- JUL-1990 

238 

1500 

28.97 

0.1245 

1 .00314 

2 

10-SEP- 1990 

574 

1500 

28.88 

0.1247 

1 .00274 

2 

1 8 - JUN- 1 990 

60 

2000 

28.80 

0.1241 

1 .00255 

2 

24- JUL- 1 990 

238 

2000 

28.89 

0.1242 

1.00240 

2 

10-SEP- 1990 

574 

2000 

28.72 

0.1243 

1 .00229 

2 

18- JUN- 1990 

60 

2500 

28.73 

0.1241 

1 .00238 

2 

24- JUL- 1990 

238 

2500 

28.77 

0.1241 

1 .00237 

2 

10-SEP- 1990 

574 

2500 

28.62 

0.1242 

1 .00232 

2 

18- JUN- 1990 

60 

3000 

28.68 

0.1243 

1 .00197 

2 

24- JUL- 1990 

238 

3000 

28.72 

0.1243 

1 .00192 

2 

10-SEP- 1990 

574 

3000 

28.56 

0.1244 

1.00200 

2 

1 8- JUN- 1990 

60 

3500 

28.91 

0.1242 

0.99983 

2 

24 - JUL- 1 990 

238 

3500 

28.93 

0.1241 

0.99983 

2 

10-SEP- 1990 

574 

3500 

28.84 

0.1243 

0.99975 



TABLE 6 

Average Sample Mapping Function Coefficients 


DETECTOR YDEF 

SO 

RMS(SO) 

B 

RMS(B) 

E 

RMS(E) 

1 

500 

30.13 

0.09 

0.1260 

0.0001 

1 .0017 

0.0001 

1 

1000 

30.65 

0.09 

0.1258 

0.0001 

1 .0020 

0.0001 

1 

1500 

31 .12 

0.07 

0.1257 

0.0001 

1 .0020 

0.0002 

1 

2000 

31.45 

0.06 

0.1256 

0.0001 

1 .0023 

0.0001 

1 

2500 

32.05 

0.05 

0.1257 

0.0001 

1 .0019 

0.0001 

1 

3000 

32.58 

0.06 

0.1255 

0.0001 

1.0023 

0.0001 

2 

1000 

28.91 

0.13 

0.1249 

0.0001 

1.0045 

0.0001 

2 

1500 

28.94 

0.05 

0.1246 

0.0001 

1.0029 

0.0002 

2 

2000 

28.80 

0.08 

0.1242 

0.0001 

1.0024 

0.0001 

2 

2500 

28.71 

0.08 

0.1241 

0.0001 

1.0024 

0.0001 

2 

3000 

28.65 

0.08 

0.1243 

0.0001 

1.0020 

0.0001 

2 

3500 

28.90 

0.05 

0.1242 

0.0001 

0.9998 

0.0001 


TABLE 7 

Prelaunch (1984) Sample Mapping Function Coefficients 


DETECTOR YDEF 

SO 

B 

E 

1 

500 

30.94 

0.1264 

1.0018 

1 

1000 

31.32 

0.1262 

1.0020 

1 

1500 

31.57 

0.1262 

1.0022 

1 

2000 

31.84 

0.1263 

1 .0020 

1 

2500 

32.25 

0.1264 

1 .0015 

1 

3000 

32.80 

0.1264 

1.0010 

2 

1000 

28.58 

0.1253 

1.0045 

2 

1500 

28.78 

0.1251 

1 .0028 

2 

2000 

28.75 

0.1249 

1.0023 

2 

2500 

28.80 

0.1248 

1.0020 

2 

3000 

28.83 

0.1249 

1 .0018 

2 

3500 

29.25 

0.1248 

0.9994 


ATTACHMENT 2 


Analysis of the GHRS Detector 
Stability Mini funct i ona Is 


D. J . Li nd 1 e r 
July 15, 1991 


ABSTRACT 

The detector stability mi n i f unc t i ona 1 was run 10 times for both GHRS 
detectors during the period of November 12, 1990 through June 9, 1991. 
Results show no changes in the focus or the geometric properties of 
the detectors over the 7 month period. The deflection step size 
remained stable to within 1 percent and the dark count remained at 
approximately 0.01 count s /s e cond/d i ode for both detectors. The only 
significant change was a decrease of observed flat field lamp count 
rate. The count rate is decreasing a rate of approximately 7 percent 
per year. It has not yet been determined whether this decrease is a 
result of a detector sensitivity loss or a decrease in the light 
output of the lamps. 


1 . INTRODUCTION 

The stability mi n i f unc t i ona 1 (proposal 2924) consists of internal 
detector observations including: 

1) A 70 second dark count observation 

2) An observation of central region of the photocathode 
illuminated by the flat field lamp. 

3) A focus check consisting of 231 observations at the 
photocathode mask edge at y-deflection 2048 and 
varying x-deflections. 

4) Left and right photocathode edge scans at y-de f 1 ec t i on 
2048 for the sample mapping function coefficient 
comput a t i on . 

5) Top and bottom photocathode mask edge scans for the 
line mapping function computation. 

2. DARK COUNT NON I TOR I NG 

The 70 second dark count observation is used to monitor any 
significant increases in the dark rate and find any diodes which 
become excessively noisy. The average dark rates in 
counts/second/diode (c/s/d) over the seven month period are shown in 
tables 1 and 2. Both detectors remained at approximately the 0.01 
c/s/d level with the detector 1 having a slightly lower level than 
detector 2. None of the observations showed any anomalous noise 
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event s or diodes. 


3 . DETECTOR FOCUS 

Tables 1 and 2 show the detector focus of the seven month period. The 
focus value tabulated is the ful 1 -width-half -maximum of a Gaussian fit 
to the differential edge profile as computed by routine FOCUS_OGMP . 
It is given in deflections units Cl deflection unit = 6.25 microns). 
The focus for both detectors has remained stable over the seven month 
period. 


4. FLAT FIELD LAMP 

The total flat field count rates in seven substep bins observing the 
central region of the photocathode are tabulated in tables 1 and 2. 
Normalized plots of the count rates versus time are shown in figures 1 
and 2. Both detectors show the count rate decreasing with time at a 
rate of approximately 7 percent per year. From this data it can not 
be deduced whether the detector sensitivity is decreasing or the lamp 
output is decreasing. The same type of analysis using spectral 
calibration lamp data show a similar decrease. Figures 3 through 7 
show the decrease in the observed count rate for the first order 
spectral calibration 1 amp mi ni f unc t i ona 1 s . Different plotting symbols 
are used to show different carrousel positions. 


5. PHOTOCATHODE MAPPING FUNCTION 

Tables 3 and 4 show the photocathode mapping functions coefficients 
computed from the detector mini func t i onal s . Results indicate that the 
deflections stepsizes (coefficients A and B) are stable to better than 
one percent. No time variations in the coefficients are evident. 


Table 1 

Detector 1 Stability 


Date 

Focus 

(def. units) 

Average 
Flat Field 

Dark Rate 
(c/s/d) 

12 -Nov -90 


6.71 

4122 

.009 

01 -Dec -90 


6.70 

4110 

.006 

15 -Dec -90 


6.62 

4085 

.008 

14- Jan-91 


6.60 

4070 

.006 

25- Jan-91 


6.66 

4063 

.005 

15-Feb-91 


6.68 

4064 

.007 

28-Feb-91 


6.56 

4052 

.013 

12 -Mar -91 


6.64 

4033 

.007 

2 8 -Ma r - 9 1 
09- Jun-91 

6.63 4011 

Table 2 

Detector 2 Stability 

.004 

.006 

Date 

Focus 

Average 

Dark Rate 


(def. 

units) 

Flat Field 

(c/s/d ) 

12 -Nov -90 


6.56 

1324 

0.010 

26 -Nov - 90 


6.57 

1325 

0.010 

10-De c - 90 


6.55 

1317 

0.009 

17- Jan-91 


6.52 

1311 

0.015 

23-Jan-91 


6.38 

1309 

0.013 

13-Feb-9l 


6.42 

1306 

0.009 

28-Feb-9l 


6.36 

1303 

0.008 

1 2 -Ma r - 9 1 


6.60 

1303 

0.007 

29-Mar -91 


6.56 

1301 

0.009 

09 - Jun - 9 1 


6.45 

1278 

0.010 



Table 3 

Detector Min i func t i ona 1 - Line Mapping Function 


ET 

date 

ID 

1 

12-NOV- 1990 

880 

1 

2 -DEC- 1990 

1125 

1 

15 -DEC- 1990 

1200 

1 

14- JAN- 1991 

1552 

1 

25 -JAN- 1991 

1701 

1 

15 -FEB- 1991 

2429 

1 

2 8 -FEB- 1991 

2844 

1 

1 2 -MAR- 1991 

2956 

1 

2 8 -MAR - 1991 

3400 

1 

9- JUN-1991 

4103 

2 

12-NOV- 1990 

891 

2 

26-NOV- 1990 

1028 

2 

10-DEC- 1990 

1158 

2 

17- JAN- 1991 

1630 

2 

23 -JAN- 1991 

1690 

2 

13-FEB-1991 

2336 

2 

28 -FEB- 1991 

2855 

2 

12 -MAR- 1991 

2945 

2 

29 -MAR- 1991 

3464 

2 

9-JUN-1991 

4114 


ROOTNAME LO 


Z0E90108T 213.98 
Z0E90208T 214.02 
Z0E90308T 214.06 
Z0E90408T 213.95 
Z0E90608T 214.04 
Z0E90508T 214.10 
Z0E90708M 214.04 
Z0E90808T 213.97 
Z0E90908T 213.99 
Z0E90B08M 213.85 

Z0E90C08T 224.59 
Z0E90D08T 224.63 
Z0E90E08T 224.60 
Z0E90F08T 224.83 
Z0E90H08T 224.68 
Z0E90G08T 224.89 
Z0E90I08M 224.84 
Z0E90J08T 224.81 
Z0E90K08T 224.84 
Z0E90M08T 224.72 


A 

YUPPER 

YLCMER 

0.1259 

348.2 

3764.0 

0.1262 

351 .9 

3759.6 

0.1261 

349.9 

3761.0 

0.1259 

349.3 

3763.4 

0.1262 

352.0 

3759.3 

0.1265 

355.9 

3754.4 

0.1262 

351.9 

3759.3 

0.1258 

347.6 

3764.9 

0.1261 

351.2 

3760.8 

0.1263 

355.0 

3759.3 

0.1237 

232.2 

3708.7 

0.1237 

231.7 

3708.5 

0.1237 

231 .8 

3708.9 

0.1237 

230.4 

3706.6 

0.1237 

232.4 

3707.2 

0.1238 

232.0 

3704.3 

0.1240 

235.4 

3701.9 

0.1237 

231.1 

3706.4 

0.1237 

231 .0 

3705.9 

0.1240 

236.4 

3702.9 


Table 4 

Detector Mi ni f unc t i ona 1 - Sample Mapping Function 


DET date ID1 ROOTNAME YDEF SO B E 


1 1 2 -NOV- 1 990 876 Z0E90104T 2048 31.51 0.1257 1.00147 
1 2 -DEC- 1990 1121 Z0E90204T 2048 31.82 0.1259 1.00074 
1 15-DEC- 1990 1196 Z0E90304T 2048 31.62 0.1257 1.00145 
1 14- JAN-1991 1548 Z0E90404T 2048 31.51 0.1256 1.00183 
1 25- JAN-1991 1697 Z0E90604T 2048 31.65 0.1258 1.00144 
1 15-FEB- 1991 2425 Z0E90504T 2048 32.08 0.1261 1.00075 
1 28 -FEB- 1991 2840 Z0E90704M 2048 31.76 0.1258 1.00135 
1 12 -MAR -1991 2952 Z0E90804T 2048 31 . 68 0.1256 1.00144 
1 2 8 -MAR- 1991 3396 Z0E90904T 2048 31.81 0.1258 1.00126 

1 9-JUN-1991 4099 Z0E90B04M 2048 32.02 0.1260 1.00100 

2 12-NOV- 1990 887 Z0E9OC04T 2048 28.65 0.1240 1.00178 
2 26-NOV- 1990 1024 Z0E90D04T 2048 28.59 0.1241 1.00187 
2 10-DEC- 1990 1154 Z0E9OEO4T 2048 28.67 0.1242 1.00184 
2 17- JAN- 1991 1626 Z0E90F04T 2048 28.70 0.1241 1.00179 
2 23- JAN- 1991 1686 Z0E90H04T 2048 28.62 0.1241 1.00183 
2 13-FEB-1991 2332 ZOE90G04T 2048 28.42 0.1240 1.00191 
2 28 -FEB- 1991 2851 Z0E90104M 2048 28.76 0.1243 1.00181 
2 1 2 -MAR- 1991 2941 Z0E90J04T 2048 28.55 0.1243 1.00173 
2 2 9 -MAR -1991 3460 Z0E90K04T 2048 28.34 0.1241 1.00185 
2 9-JUN-1991 4110 Z0E90M04T 2048 28.66 0.1244 1.00160 



DETECTOR 1 FLAT FIELD LAMP 



6)DJ )UnOO pSZjlDCUJON 


48200 48220 48240 48260 48280 48300 48320 48340 48360 

Modified Julion Dote 



DETECTOR X FLAT FIELD LAMP 



©}dj junoo P9 Z»|dujjon 


Modified Julion Dote 


SC2 



48050 48100 48150 48200 48250 48300 48350 





48100 48150 48200 48250 48500 48550 48400 48450 






48100 48150 48200 48250 48300 48350 48400 48450 





xnu dszudlajjoki 


48100 48150 48200 48250 48300 48350 48400 48450 





ATTACHMENT 3 



GHRS Geomagnetic Image Motion 


Don J . Li nd 1 e r 
July 31 , 1991 


ABSTRACT 

We have me asured the GHR S i ma g e mo tion do to the Earth s ma g n e t i c 
field and found that the motion at the diode array is approximately 50 
microns (one diode)/Gauss foT detector 2 and approximately 10 microns 
(0.2 diodes ) /Gauss for detector 1. This motion can cause errors in 
target acquisition, errors in measured wave 1 eng ths , and can result in 
ioss of spectral resolution. The impact of the magnetic motion can be 
minimized bv limiting spectral observations to less than 5 minutes and 
bv using the Double Locate mode of target acquisition. 

1.0 INTRODUCTION 


The GHRS OY'SY program has no formal tests for geomagnetic image 
motion. The best data for estimating the motion are the series of 
deflections calibrations (DEFCALs) taken during the coarse and fine 
a 1 i g nme nt tests including; 


Proposal Date Number of DEFCALs 


2988 

13- SEP- 1990 

12 

2765 

3- JAN-1991 

15 

2765 

14-FEB- 1991 

15 


Each test contained DEFCALs distributed over a 6 to 7 hour period. 
Each DEFCAL observation consists of an onboard measurement of the 
locations of the spectral calibration 1 amp apertures observed with 
mi r r o r N1 . Th e GHRS flight sof twa re c ompu tes the differences of the 
aperture locations from their nominal location to the nearest 
deflection unit ( 6 . 25 mi c rons ) . These d i f f e r ence s we re used to track 
the X and Y image motion over the course of each test. 


A second set of observations that were used to estimate the 
geomagnetic image motion were two observations designed to measure the 
carrousel stability over a period of one hour (proposal 2168, 
1 5-Oc t - 1990 ) . The data consisted of 150 spectral lamp observations 
taken over an hour in each Echelle mode. The carrousel position 
remained fixed during each sequence. We measured the motion in the 
x-direction (dispersion direction) of the digicons by cross 
correlation of each of the spectTa with the first spectr um of the 
sequence . 
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To separate magnetic motion from other sources of image motion (e.g. 
thermal), we fit our results to the following model; 

Motion =F(t)+S*B 

where F( t ) (a quadratic function of time, t) is used to model 
non-geomagnetic motion, B is the ear ths magnet i c field component in 
the appropriate direction, and S is a factor giving the geomagnetic 
motion in microns per Gauss. The coefficients of the quadratic 

function and S are solved simultaneously by a least squares fit. 

2.0 Results 

The results for the coarse and fine alignment DEFCAL sequences are 
shown in figures 1 through 6. The diamonds show the average 
difference of the two spectral calibration lamp apertures from their 
nominal positions. The dotted lines show the quadratic functions used 
to model non-geomagnetic motion. The solid lines show' the models 
including both the geomagnetic and non-geomagnetic motion. The 
c ompu ted £ e oma gnetic sensitivity factors we r e , 

Date Direction Detector S (mi c rons /gaus s ) 


13- SEP- 1990 X 

1 3- SEP- 1990 V 
03- JAN- 1991 X 
03- JAN- 1991 Y 

14 - FEB- 1991 X 
14 -FEB- 1991 Y 

Expect for the first measurement, the observed magnetic motion for 
detector 2 is approx ima te ly one diode per Gauss. The quality of the 
present data makes it difficult to conclude that the magnetic 
sensitivity is variable. However, tests run on the Faint Object 
Spectrograph confirm that its geomagnetic sensitivity varies. 

Figure 7 show's the spectral motion in the x-direction for data taken 
w'ith Echelle B (detector 2) over a period of one hour. The 
sensitivity factor (39 mi c rons /Gaus s ) for this observation is 
consistent" with the DEFCAL observations. Figure 8, which show's the 
results for echelle A, indicates that detector 1 does not have a 
significant sensitivity to the Earth’s magnetic field. 



3.0 Reconmenda t i ons 

We make the following reconmenda t i ons 
geomagnetic image motion. 


to combat the problems with 
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1) Limit detector 2 spectral observations to 5 minutes or less. There 

is no readout noise with digicons and the readout overhead time 
is minimal. Multiple readouts can be registered and coadded on 
the ground. 

2) Use the double locate option in target acquisition. The success of 

vout acquisition will depend on the mi n imi za t i on of the time 
between the deflection calibration and the target locate or 
peakup. This is particularly important if you intend to place 
the target in the small science aperture. 
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ATTACHMENT 4 


GHRS OV/SY Wave 1 eng th Calibration 


Don J . L i nd 1 e r 
14 Augus t 1991 


ABSTRACT 

We have analyzed the GHRS OV/SY spectral calibration lamp observations 
and found that an average d i spe r s i on coef f i c i en t table with linear 
thermal motion models can be used to compute wave 1 eng ths for science 
observations to better than one photocathode sample unit (50 microns). 
If spectra] calibration lamp observations are taken along with the 

science observations, the errors can be decreased to .17 sample units 

(one sigma error), which equals 0.5 km/ s in the echelle modes. 

1.0 I nt roduc t i on 

We used the following process to construct the dispersion coefficient 
table required by the standard GHRS reduction routine, CALHRS. 

1) Compute the dispersion relation for multiple carrousel 

positions and thermal conditions using routine WAVECAL 

(section 2.0). 

2) Derive each dispersion coefficient as a polynomial function of 
carrousel position (section 3.0). 

5) Construct a thermally induced motion model by fitting the 

residuals of the fit in step 2 as a linear function of a 
selected GHRS thermistor reading (section 4.0). 

8) Generate a dispersion coefficient table on a finely spaced 
carrousel position grid (section 5.0). 

The above procedure was repeated for each grating mode. 


2.0 The GHRS dispersion relation. 

The dispersion equation used by the GHRS standard reduction routine 
CALHR S is given by : 


s 


= a 0 -i- a 1 *m*w + a 2 4- a 2 *m + r 

a5*m*m*\v + a6*m*w*w + aT^fm’mWw 


s - is the photocathode sample position 
w - is the wa v e 1 e n g t h 

m ^ is the spectral order (l for first order gratings) 
aO, al, ... a7 - are dispersion coefficients which vary 
with carrousel position 


We computed the dispersion coefficients using observations of the 
internal spectral calibration lamps SCI and SC2 with GHRS/IDL routine 
WAVECAL. WAVECAL finds the positions of spectral lines in the a 
pla t inum^neon hollow cathode lamp line list (Reader, 1990). WAVECAL 
then computes the dispersion coefficients by a least squares fit which 
minimizes the sum of the squares of the difference between the 
observed spectral line positions and .those .computed by the dispersion 
relation. 


We have found that the following non-zero coefficients are required 
for each of the gratings: 


G- 1 

(G140M) 

aO , 

a 1 , 

a 2 


G-2 

(G160M) 

aO, 

al , 

a2 , 

a7 

G-3 

(G200M) 

aO, 

al , 

a2, 

a 7 

G-4 

(G270M) 

aO, 

al , 

a2 , 

a 7 

G-5 

(G140L) 

aO , 

al , 

a2 


E-A 

(ECH-A) 

aO. 

a 1 , 

a 2 , 

a4 

E-B 

(ECH-B) 

aO, 

al , 

a 2 , 

a4 , a 7 


The cubic coefficient, a7, is required only for the detector 2 
gratings and coefficient a4 is redundant with the al coefficient for 
the first order gratings where the spectral order, m, equals 1. 
Analysis of the residuals of the dispersion coefficient fits indicate 
that coefficients a3 and a5 are not required. 

Experiments have shown that the cubic coefficient, a7, can not be 
accurately fit for every carrousel position. Th e spectral calibration 
lamp spectra have an inadequate number of spectral lines or poorly 
spatially distributed spectral lines at many of the carrousel 
positions. The cubic term does not, however, vary significantly with 
carrousel position. This allows us to compute a constant a7 for each 
grating by fitting all carrousel positions simultaneously. For 
example. Figure 1, shows the residuals (observed minus fitted spectral 
line positions) for all G-2 observations taken at multiple carrousel 
positions when no cubic term is used. A cubic polynomial was then fit 
to these residuals to give an a7 coefficient applicable to all 
carrousel positions for the grating. We then r e r an WAVECAL on each 
observation with a fit to the non-cubic coefficients while holding the 
cubic term, a7, fixed at the the value computed from all carrousel 



positions. Figure 2 shows the same plot as figure 1 when cubic term 
is j nc 1 ud ed . 


3.0 Dispersion coefficient variations with carrousel position. 

The dispersion coefficients computed by WAVECAL do not vary smoothly 
with carrousel position (Cushman, Ebbets, and Holmes, 1986). This is 
particularly true of the aO and al terms. These unpredictable 
variations between carrousel positions make it difficult to analyze 
thermal variations and to compute dispersion coefficients for 
carrousel positions without a corresponding spectral calibration lamp 
observation. The solution to this problem is to transform the 
coefficients to a different coordinate system in which they vary 
smoothly. The following coordinate system satisfies this requirement. 

s = fO + fl*x + f2*x*x + f3*x*x*x + f'4*y 

x = m*w - me -wc 
y = w - wc 

me is the central order (42 for E-A, 25 for E-B, and 
1 for the first order gratings) 
wc is the central wavelength at the given carrousel position 
computed by the carrousel equation: 

A C - R 

wc = * S I N ( - -) 

me 10430.378 

R is the carrousel position 

A and C are coefficients which vary with grating mode and 
are determined by a non-linear least squares fit 
to the observed wavelengths at photocathode sample 
position 280.0. 

The coefficients fO, . . . , f 1 are calculated .from the dispersion 
coefficients aO al at carrousel position R by: 

k = me *wc 

fO = aO + al*k + a2*k*k + a7*k*k*k + a4*wc 
fl - al - 2*a2*k + 3*a7*k*k 
f2 « a2 + 3 * a 7 * k 
f 3 = a 7 
f 4 = a 4 

We can now fully characterize the dispersion relation for all 
carrousel positions of a grating by fitting polynomials of carrousel 
position to the f 0 , fl, f2, f3, and f4 terms: 


fO = F00 - F01 *R - F02*R*R 
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fl = FI 0 + FI 1 *R + FI 2 *R*R 


f 4 = F40 - F4 1 *R 4- F42*R*R 

wtere R is the carrousel position and Fi j are t he coe f f i c i ent s of the 
polynomials. The derived values of Fij for all grating modes are 
given in table 1. 


4 . 0 Spe c t r a 1 Motion 

The residuals of the polynomial fits to the fO, fl, ... and f4 terms 
in section 3.0 can be used to analyze changes in the dispersion 
relation with respect to temperature, time and the earths magnetic 
field. Since we had a large number of wavelength calibration 
observations for each grating mode., we treated the values of the 
coefficients, Fij, representative for average conditions. The 
difference of the dispersion coefficients for each individual spectral 
calibration lamp observation from the values of the polynomial fits 
(which we will call the predicted values) can be attributed to thermal 
motion, geomagne t i ca 1 ly induced motion, carrousel repeatability, 
instrument aging, ect. The largest differences were found in the fO, 
or constant, term. Changes in this term represent movement of the 
entire spectral format. We have correlated the offsets of the fO 
coefficient with temperature, time, and the relevant component of the 
Earth’s magnetic field vector. The largest source of variation in the 
fO term is caused by thermal effects. Table 2 gives the results of 
linear fits to the offsets for 22 different temperatures listed in 
table 3. The RMS of the offsets of fO from the predicted values and 
the maximum of the absolute value of the offsets are given in the row 
titled ”No Fit". The other rows give the RMS and maximum of absolute 
value of the residuals of a linear fit of the offsets to the specified 
thermistor values. The RMS is an estimate of the correlation with the 
specified thermistor reading. . If the. RMS is near that of the "No Fit" 
value, no significant correlation exists. The results of the best 
correlations are shown in Figures 3 through 9. The offsets of fO from 
the predicted values for each spectral calibration observation are 
shown as plus marks (4). The solid lines show the linear thermal 
mode 1 s . 

The values of the RMS columns in table 2 are useful for determining 
the accuracy of the ’wavelengths assigned by CALHRS . If an average 
dispersion coefficient table is used without application of a thermal 
model, the "No Fit" rows in table 2 give an indication of the typical 
errors. For example, the typical error for grating 4 would be 2 
sample units (100 microns) and in some cases could be almost 4 sample 
units ( 200 microns). Using the linear thermal model option of CALHRS 
reduces the typical grating 4 error to 0.25 sample units. 


Table 2 

Linear Thermal Model for Spectral Motion in the Dispersion direction 
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Table 3 


Table 2 Thermistor Locations 


ZR1LTA 

ZRILTB 

ZDETT1 

ZDETT2 

ZDEBTF 

ZDEBTR 

ZPABT1 

ZPABT2 

ZMEBTl 

ZMEBT2 

ZF1AT 

ZFIBT 

ZFICT 

ZCST 

ZSCT1 

ZSCT2 

ZHVPST1 

ZHVPST2 

ZDT11 

ZDT12 

ZDRT 

ZOBBT 


RIU A temperature 

RIU B temperature 

Detector 1 temperature 

Detector 2 temperature 

DEB front post -amp temperature 

DEB rear post -amp temperature 

Detector 1 preamp assembly box temperature 

Detector 2 preamp assembly box temperature 

.Main electronic box 1 temperature 

Main electronic box 2 temperature 

Fixture interface A temperature 

Fixture interface B temperature 

Fixture interface C temperature 

Carrousel stator temperature 

Spectral calibration lamp 1 temperature 

Spectral calibration lamp 2 temperature 

Detector 1 high voltage power supply temperature 

Detector 2 high voltage power supply temperature 

Detector 1 shield temperature 

Detector 2 shield temperature 

Digicon radiator temperature 

Optical bench bulkhead temperature 


We next investigated the variability of the fO term over time. Only 
grating 5 showed any correlation with time. Figure 10 shows the 
offsets (with motion from our thermal motion model removed) versus 
Modified Julian Date for grating 5. The correlation is not convincing 
but it does indicate that time variations should be monitored in the 
future . 

A final source of motion we investigated is the detector 2 sensitivity 
to variations in the Earth’s magnetic field (Lindler, 1 9 9 1 ) . Figure 
11 shows a plot of the fO offsets versus the geomagnetic field 
component in the direction of dispersion. Some correlation is 
evident. The slope of the least squares linear fit (shown as the 
solid line) is approximately 0.75 sample units per Gauss. This is 
consistent with other measurements of the geomagnetic sensitivity. 

f 1 also showed significant variations.. Change s .in this term represent 
a plate scale change. Figures 12 through 18, show the differences in 
the linear dispersion term from the predicted value as a percent 
change. The changes are plotted against the thermistor reading which 
gave the best correlation. A change in dispersion of 0.05 percent can 
result in a 0.25 sample unit error over the 500 diode array 
(500*0.0005). CALHRS presently has no provision for handling a 
thermal model of changes in the linear dispersion term. Users with 
observations reduced with an average dispersion coefficient table 
should be aware that not only is there a zero point shift in their 
wavelengths but also changes in dispersion. 


5.0 Dispersion coefficients for arbituary carrousel positions. 

The present implementation of CALHRS uses a table of dispersion 
coefficients tabulated on an arbitrary grid of carrousel positions. 
If wavelengths are required for an untabulated carrousel position, 
CALHRS linearly interpolates the coefficients between the two closest 
carrousel positions. With the s i gni f i can t. improvement s to the thermal 
model. this linear interpolation becomes one of the most significant 
remaining errors. Using the polynomial model for the dispersion 
coefficients from section 3.0, this error can be made negligible by 
creating a dispersion coefficient table on an arbitrarily small 
carrousel grid. We chose a grid with dispersion coefficients 
tabulated at every carrousel position divisible by 4. This 
effectively eliminated interpolation errors at the cost of 
substantially increasing the size of the dispersion coefficient table 
and thus increasing the CALHRS execution time. 



6.0 CALHRS reference files 


The final results of this report are tabulated in two tables for use 
by CALHRS. Table DC__004.TAB contains the dispersion coefficients for 
all grating modes, tabulated at carrousel positions divisible by 4. 
The coefficients were generated using the polynomial coefficients 
listed in table 2. A thermal coefficient table, T\L_002.TAB, was 
generated by selecting the thermistor readings for each grating mode 
which gives the best correlation with spectral motion. The thermal 
motion coefficients in TM_002 are listed in table 3. 

Table 3 

Linear thermal motion coefficients 

Grating Thermistor Motion 

( s amp les/degree ) 


G-l 

ZDEBTR 

-0.03 

G-2 

ZRIUTA 

-0.09 

G-3 

ZFIAT 

-0.45 

G-4 

ZCST 

-0.69 

G- 5 

ZFIBT 

0.26 

E-A 

ZCST 

-0.37 

E-B 

ZDEBTR 

-0.43 


7.0 Errors in the wavelengths assigned by CALHRS. 

There are a large number of sources which produce errors in the 
wavelengths generated by CALHRS. Many of the these errors can be 
significantly reduced by using proper observing techniques or careful 
reduction techniques. All errors described in this section are quoted 
in photocathode sample units. One sample unit equals 50 microns which 
is the separation between neighboring diodes. 

7 . 1 Image motion. 

The largest sources of wavelength error are thermal motion and 
carrousel repeatability. Table 4 gives the RMS error and the maximum 
error observed to date for each of the grating modes. The errors 
occur when performing reductions using the average dispersion 
coefficient table (DC_004) and no thermal motion correction. 



Table 4 

CALHR S errors using the average 
dispersion coefficient table (DC_004) 


Grating 

RMS error 

Maximum error 


( Samp 1 e units) 

( s amp 1 e units) 

G-l 

.4 

1 .0 

G - 2 

.5 

1 .3 

G - 3 

1 .6 

3.3 

G-4 

2.0 

3.7 

G- 5 

.9 

2.0 

E-A 

.9 

1 .9 

E-B 

.5 

1 .0 

The errors can be reduced to those shown. in table 5 by apply 

simple linear thermal 

motion mode 1 . de.s c r i-bed by TM_002 . 


Table 5 


CA1HRS errors us 

ing the average 

dispersion coefficient 

table (DCL004) and thermal model 

table CIM_002) 

Grating 

RMS error 

Maximum error 


( Sampl e units) 

( s amp 1 e units) 

G-l 

.2 

.4 

G-2 

.3 

.7 

G-3 

.4 

.9 

G-4 

.2 

.7 

G-S 

. S 

1.7 

E-A 

.2 

.4 

E-B 

.3 

.6 


Thermal effects and other sources of motion errors can be further 
reduced by observing the spectra] calibration lamp (WWEC4L) before 
and/or after the observation of the target. The thermal motion and 
carrousel repeatability errors are negligible if the WAVECAL is taken 
at the same carrousel position as the science observation. The only 
concern is for the short-term motion between the time of the spectral 
lamp and science observation, typically about 0.2 to 0.4 sample 
units /hour . 

7.2 Dispersion changes. 

Unlike thermally induced motion, changes in dispersion give errors in 
differences of wave 1 eng ths within the same observation. The linear 
dispersion at a carrousel position can vary by 0.1 percent. This can 
result in a 0.5 sample unit error when wavelengths measured at the two 
ends of the diode array are compared. This error can be effectively 



eliminated by using the dispersion coefficients from a WWEC4L taken 
before or after a science observation. 

7.3 Aperture Offsets. 

The wavelength offsets between the spectral calibration lamp apertures 
and the small science aperture (SSA) were cal ibra ted us ing prelaunch 
data taken during phase 6 (1984) calibration at Ba 1 1 Aerospace. The 
estimated (one sigma) error in this calibration is 0.1 sample units. 
There is no way for an observer to decrease this error. However, it 
is important to note that the error repeats. If you observe with the 
same carrousel position at two different times, the aperture offset 
error is the same. A difference of two measured wavelengths from one 
observation to the other will not show an aperture offset error. 

The calibration between the SSA and the 1 a r ge . s c i ence aperture (LSA) 
has not yet been completed. CALHRS presently computes SSA wave 1 eng t hs 
for objects in the LSA. Errors resulting from this lack of 
calibration are probably on the order of 1 to 2 sample units. 

7.4 Target centering. 

CALHRS computes wave 1 eng ths assuming that an object is centered in the 
target aperture. If an object is not centered, wavelength errors will 
occur. Mi scenter ing can result from errors in the deflection 
calibration, target locate or peakup, thermal or geomagne t i ca 1 ly 
induced image motion between the defcal and peakup, errors in the GHRS 
onboard values for the aperture locations, errors in the slew between 
the LSA and SSA, and spacecraft drift/jitter. Because of the size of 
the SSA and the size of the spherically abberated stellar image the 
one sigma error f r om mi s cent e r i ng in the SSA is on the order of 1/4 of 
a sample unit. The best measurement of miscentering of a target in 
the SSA is the comparison of the flux measurement between the LSA and 
SSA: if the SSA flux is approximately 1/4 of the LSA flux, the object 
is well centered. The centering error should be significantly reduced 
once the new on-board SSA flux peakup routine is implemented 

7.5 Geomagnetic image motion. 

The detector 2 sensitivity to the Earth’s magnetic field can result in 
loss of spectral resolution and wavelengths errors up to 0.5 sample 
units. To avoid loss of spectral resolution, detector 2 spectral 
observations should be limited to short exposures (i.e. 5 minutes or 
less). If it is found that the detectors sensitivity to the Earth’s 
magnetic field does not vary, it will be easy to model the geomagnetic 
motion and remove it. If the sensitivity does vary (as shown to be 
the case f ot the Faint Object Spectrograph), a WAVEC^L taken before 
and after the science observation could be used to calibrate the 
sensitivity value at the time of the science observation. 


7.6 The dispersion model. 


Dispersion coefficients are computed by least squares fits of the 
dispersion model to spectral line positions measured in spectral 
calibration lamp observations. A measure of the errors in the 

dispersion coefficients is the RMS of the differences between the 
observed spectral line positions and the positions computed by the 
dispersion model. We found typical RMS errors of 0.04-0.08 sample 
units for the first order gratings and 0.1-0.15 sample units for the 
echelle modes. The RMS errors in the echelle modes are slightly 
larger because of the inadequacies in the dispersion model when the 
entire format (multiple spectral orders) was fit s imu 1 a t aneous 1 y . The 

RMS residuals for the echelle modes can be reduced to less than 0.1 

wti en single orders are fit. 

8.0 Ultimate wavelength precision. 

If proper observing techniques are used, what is the ultimate 
wavelength precision of the GHRS? I f WAVECALs are taken before and 
after the science observation in the SSA, errors from image motion, 
carrousel repeatability, and geomagnetic motion can be rendered 
negligible. The remaining errors are: 

1) Errors in the dispersion coefficients (0.1 sample units) 

2) Errors in the wavelength offsets between the spectral 

calibration lamp and small science aperture (0.1 sample units) 

3) Errors caused by improper centering of the object in the SSA 

With the present state of instrument calibration, errors in the 

centering can result in wavelength errors of approximately 0.25 sample 
units (one sigma). Once the SSA flux peakup option becomes available, 
this error will probably be reduced to approximately 0.1 sample units 
(one sigma). Combining these three sources of errors (0.1, 0.1, and 

0.1) gives the residual error of 0.17 (the square root of [3 x 0.1 
squared]) which equals approximately 0.5 kin's in the echelle modes. 
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Figure Capt i ons 


figure 1: Wavelength calibration residuals of the observed grating 2 

spectral line positions from those of least squares quadratic 
dispersion model. Residuals are in 50 micron sample units. The 
plot shows the need for a cubic dispersion term. 

figure 2: Same as figure 1 except that a cubic dispersion model was 

used . 

figures 3-9: Linear thermal motion models for each grating. The plus 

marks (+) show the offset of each spectral lamp observation from 
an average dispersion model for the grating plotted against the 
GHRS thermistor value with the best correlation. The solid lines 
are linear least squares fits. Offsets are in 50 micron sample 
units and temperatures are in degrees Celsius. 

figure 10: Grating 5 spectral motion with time. The plus marks (+) 

show the offset of each spectral lamp observation from an average 
dispersion model for the grating plotted against time. The solid 
line is a linear least squares fits. Offsets are in 50 micron 
s amp 1 e units. 

figure 11: Detector 2 spectral motion plotted against the conponent 

of the Earths magnetic field in the dispersion direction. The 
plus marks (+) show the offset of each grating 2, 3 or 4 spectral 
lamp observation from an average dispersion model for the 
grating. The solid line la a linear least squares fits. Offsets 
are in 50 micron sample units. 

figures 12-18: Changes in the linear dispersion with temperature. 

The plus marks (+) show the change of the dispersion for each 
individual spectral calibration lamp observation from an average 
dispersion model for the grating. The solid line shows a the 
least sqaures fit of the changes with the GHRS thermistor giving 
the best correlation with temperature. The temperatures are 
plotted in degrees Celsius. 


Residuals (somple units) Residuals (sample units) 


Figure 1 
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Figure 2 
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GHR S Wa ve length Of f s e t s 
Sma J 1 Science 


Between the Large and 
Ape r t u r e s . 


Don Lmdler 
October 30, 1991 


ABSTRACT 

We have made measurements of the offsets between the large (LSA) and 
small (SSA) science apertures for GHRS gratings 2, 3, and 4. In all 
cases the measured offsets were less than one diode and the average 
offset for each grating mode was less than 1/2 diode. There w'as 
insufficient data to quantify changes in the offsets versus carrousel 
position and photocathode sample positions. . However, grating 2 
results indicate that the offset does vary with carrousel position. 


1) INTRODUCTION 


The standard GHRS reduction routine, CALHRS , has the provision for 
applying a wavelength offset to the wavelengths computed from a table 
dispersion coe f f i c i en t s . The dispersion coefficients are used to 
compute wavelengths applicable to data taken in the small science 
aperture. The wave length offsets for the other apertures are computed 
by : 

EXV = (A + B * s ) /m (l) 

whe r e , 


EXV 7 is the offset added to the sma 1 1 science aperture 
wavelength values. 

s is the photocathode sample position of the data point, 
m is the spectral order (l for the first order gratings) 


A and B are offset coefficients that are tabulated by grating 
mode, aperture, spectral order, and carrousel 
position. 


Pre launch measurement s of the offsets between the small science 
aperture and the spectral calibration lamp apertures (SCI and SC2) 
show va r i a t i ons in A and B for changes in both carrousel position and 
spectral orders. Results also show that B is non-zero. This means 
that not only is there a wavelength offset between entrance apertures 
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but there is also a plate scale chance. 

Measurements of the offsets between the SSA and the spectral 
calibration lamp apertures we r e ma de in 1984 by ill um inating the S S A 
with a platinum neon spectral calibration lamp. No measurements of 
the oflsets between the SSA and LSAwere made during the prelaunch 
cal i b r a t ion. 


Data for measurement of the offsets between the SSA and LSAwere taken 
on June 15, 1991 (proposal 2097). The data consisted of observations 
of HD93521 taken in both the large and sma 1 1 science apertures. Prior 
to each observation, wavelength calibration observations were taken so 
that any effects of thermal motion, geomagnetic induced motion, and 
carrousel repeatability could be removed. Five carrousel positions 
were observed for grating 2 and three carrousel positions were 
observed for gratings 3 and 4. 

2) DATA REDUCTION 

Fo r each carrousel po sition the foil owi ng sequence of observations wa s 
made . 


A) Spectrum Y Balance observation 

B) Waveca] observation using lamp SCI 

C) HD93521 observation through the LSA 

D) Spectrum Y- Balance observation 

E) Wavecal observation using lamp SCI 

F ) HD9 3521 observation through the SSA 

Observations B and E were used to remove any of the wave length offsets 
between observations C and F which resulted from thermal and 
geomagnetic image motion or from the carrousel motion needed for the 
spectrum Y -balance observations. The offsets between observations B 
and E and observations C and F were computed using a normalized mean 
and variance correlation. Quadratic refinement of the correlation 
values was used to determine the offset to sub-pixel accuracy. 
Observations C and D were smoothed to the same resolution before 
correlation by convolving the LSA observation with the SSA line spread 
function and the SSA observation with the LSA line spread functions. 
To minimize the effect of fixed pattern noise, the correlations of 
HD93521 were limited to regions with large spectral features. 

The offset between the SSA and LSA was then computed as the offset 
between observation C and F minus the offset between observations B 
and E . 
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3) ANALYSIS 

The offsets that we computed between the SSA and LSA are tabulated in 
iah]e 1 versus grating mode. carrousel position, and pho t oc a t hod e 
sample position. The photocathode sample position column gives the 
position of the center of the HD93521 spectral 1 eature that was 
correlated. 

The results for grating 2 show some evidence that the offsets vary 
with carrousel position. The grating 2 observations at carrousel 
position 51508 also show some evidence that the offset also varies 
with sample position. Variations with sample position would indicate 
that the plate scale or dispersion is different between the two 
apertures . 

T o i mp 1 erne n t the LSA/ SSA offset correction in CALHR S . we have elected 
to only use an average offset for each grating. The minimal amount of 
data does not allow us to construct accurate models of the variations 
of the offset with carrousel position and sample position. The 
resulting coefficients (Equation l) for use by CALHR S are then: 


Grating 

A 

B 

G 2 

-0.022 

0.0 

G 3 

-0.037 

0.0 

G - 4 

-0.036 

0 . 0 


It should be noted that there ma y be errors in these coefficients 
resulting from the target acquisition of HD93521 in the LSA and SSA. 
These errors include; statistical errors in the on-board deflection 
calibration and LSA locate, errors in the calibrated position of the 
LSA with respect to the calibration lamp apertures, and errors in the 
slew between the LSA and SSA. No peakup was performed to center the 
object in the SSA. 
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tabu : i 


Wavelength Offsets Between the 
1 . a r g e and Sma 31 Science Apertures 


Oral i n g 

Carrousel 

Position 

Photocathode 
Samp 1 e Po s i t on 

Offset 
(diodes ) 

()l f se t 
■ An s t oms 

G - 2 

49936 

244 . 

0.82 

0.058 

G- 2 

50680 

229. 

0.48 

0.034 

G-2 

51096 

164 . 

0.70 

0.050 

G - 2 

51096 

292 . 

0.44 

0.03 3 

G-2 

51 304 

417 . 

0.55 

0.039 

G-2 

51508 

85 . 

-0.08 

-0.005 

G-2 

51508 

160. 

-0.09 

-0.006 

G-2 

51508 

273. 

0.04 

0.003 

G-2 

51508 

323 . 

0.17 

0.012 

G-2 

51508 

404 . 

0.09 

0.007 

G - 3 

25056 

335. 

0.38 

0.030 

G-3 

25876 

354. 

0.59 

0.046 

G - 3 

26288 

276 . 

0.31 

0.024 

G-3 

26288 

367 . 

0.63 

0.049 

G-4 

9360 

242 . 

0.46 

0.042 

G - 4 

9 360 

317. 

0.44 

0.040 

G-4 

9832 

235. 

0.25 

0.022 

G-4 

9832 

339. 

0.41 

0.037 

G-4 

9832 

395. 

0.35 

0.031 

G-4 

9832 

483. 

0.44 

0.040 

G-4 

10256 

223. 

0.48 

0.043 

G-4 

10256 

311 . 

0.39 

0.035 
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BLOCK ITERATIVE RESTORATION OF ASTRONOMICAL IMAGES FROM 

THE HUBBLE SPACE TELESCOPE 

Don J. Lindler 

Advanced Computer Concepts 
Potomac, Maryland 


1. INTRODUCTION 

The discrete model of Unear image degradation is specified by the equation: 

b = Hx -f n (1) 

where b and x are the pixel values of the degraded and original undegraded images stacked into 
column vectors, H is a matrix constructed from the impulse response (or point spread function) 
of the degradation, and n is an unknown additive noise vector. The object of restoration is 
to determine x, given b and possibly information on the properties of n. If the point spread 
function used to construct H is not known for the given optical- detector configuration, it must 
be estimated from the blurred image, b. The point spread function is most easily estimated 

from point sources (t.e. stars) on the blurred image. 

Since H may be ill-conditioned or singular, and only the statistical properties of the noise are 
known, there are many solutions for x which satisfy equation (1). The success of a restoration 
therefore depends on the abiUty to model and apply to the restoration, known or assume 
properties of the desired solution, such as positivity or smoothness. 

Some advantages of algebraic image restoration are: 

1) The point spread function may be spatiaUy variant; 

2) If a constrained least squares method is used, the applied constraints may be varied from 
pixel to pixel to make maximum use of the known image properties; 

3) Missing or bad pixel values in the blurred images can be easily handled without attempting 
to repair their values; 

4) Noise properties can vary from pixel to pixel. 

The main disadvantage of algebraic image restoration is the size of the linear system. For a 
500 x 500 pixel image, H is a 250,000 x 250,000 matrix. Even with the most powerful computers 
available, a direct solution of the system would be impossible. In the next section, we describe 
a technique - the block iterative method, of solving large linear systems. 

2. THE BLOCK ITERATIVE RESTORATION ALGORITHM 

2 A Block Jacobi Iteration 

In most astronomical images, the point spread function has a much smaller spatial extent 
than the image, so it is appropriate to work on the image locally. We therefore divide the 
image into blocks and restore each block separately, using values from the previous iteration as 
estimates of the unblurred image values outside the block. In most instances the blurred image 


is a good choice for the starting or zeroth iteration. This type of iteration is called block Jacobi 
or group Jacobi iteration (Young 1971) and can be formulated in matrix notation as follows. 
Consider the blurred image, b, divided into m blocks of equal size 2?,, 1 = l,m: 

B 2 ^ 

D Bi - 1 B t B, +1 

n — 

^ Bm- 1 B m 

Stack the elements of each block and place them into a vector: 

B = (B l ,B 2 .---,B m ) T 

Ignoring the noise for now, we write the system as: 

HX = B 


where H is partitioned into blocks 


H „ 

H u 

• H lrn 

\ 

H 21 

H 22 

■■■ H 2m 


H m ! 

H m 2 

* ' * H mm 

/ 


and X contains the restored values, blocked in the same manner as B. If the image were divided 
into blocks of n pixels each, then the block H i: would have size n x n. The block Jacobi method 
can now’ be written as: m 

Hi.x ;* 1 = Bi- £ h„x; m 

i = 1, - • • , m, and where XJ is the stacked values for iteration r of block j. If we define the vector 
on the right hand side of equation (2) as BMOD, (i.e., the blurred image less contributions 
from outside the block as estimated from the previous iteration), the linear system for block i 
can now be written as: 

HnX r+l = BMODi (3) 

Using the block Jacobi method, we can reduce the problem to solving 

Hx — b ( 4 ) 

where H is Hu for block i; x is X t r+1 for block t and iteration r; and b is BMOD, for block t. 

The solution for block i now requires the solution of an n x n linear system. For example, 
to restore a 100 x 100 pixel image divided into m=100 blocks, each of size n x n = 10 X 10, the 
largest system to be solved would have Hu of size 100x100. Since solutions of linear systems 
require on the order of n 3 operations, the block approach compares favorably to the direct 
solution of the 10,000 x 10,000 system. For a spatially invariant point spread function, the 
problem is further reduced because Hu will be identical for all t = 1, • • - , m. 



If a constrained least squares approach is used to solve the linear system, the solution will 
converge to acceptable results even with a block size as small as the full-width-at- half-maximum 
(FWHM) of the point spread function. Overlapping the blocks (accepting only the central por- 
tion for the next iteration) can be used to speed convergence. 

2,2 Image Constraints 

The block Jacobi method reduces the restoration to solution of many smaller linear systems, 
but it does not address the ill-conditioned nature of H or the presence of noise in the blurred 
image. An ill- conditioned matrix means small changes in b, caused by noise, yield large changes 
in the solution x = H~ 1 b. In this section, we show how constrained solutions can handle these 

problems. 

In most images, the data vary smoothly except at isolated points or edges. For example, an 
image of a star field will vary smoothly, except at locations of individual stars. We can make 
use of this image properly, smoothness, by applying a constrained least squares fit. Specifically, 
we minimize a linear operator \\Qz\\ (i.e. the sum of the squares in Qx), where Q is a matrix 
designed to control smoothness or other characteristics of the image (Twomey 1963, Philips 
1962). For example, we can control smoothness in the one dimensional case by minimizing the 
second difference in the solution subject to some other constraint. If the statistical properties of 
the noise are known, we could minimize the second difference such that the norm of |iffx-b|| = 
n: that is to say, the difference of the blurred image and the solution reconvolved with the point 
spread function should have the same properties as the noise. In this case (minimize the second 
difference), Q would have the form: 



We use the method of Lagrangian multipliers, sometimes called the method of undetermined 
multipliers, to compute a solution, x, given by (Andrews 1977): 

x = ( H t H + iQ t Q)~ 1 H t b (5) 

7 is the reciprocal Lagrangian multiplier which can be selected to control the smoothness of the 
solution. Solutions using LagTangian multipliers place no restrictions on the form of Q. This 
flexibility allows the development of a variety of constraints depending on the known properties 

of the image. _ . — 

Figure 1 shows the application of this constrained least squares filter for a test case (a point 
source) with different values of 72 - The subscript 2 is used to indicate that the constraint is 
the minimum second difference. Note in figure l.c, with the largest value of 72 , noise in the 
solution has been supressed. However, the width of the point-source profile is almost as wide 
as the blurred profile. Also, some ringing in the restored profile is evident. Restored values on 
each side of the profile drop significantly below the background level. These problems result 


because the second difference is large at the location of a point source. We therefore remove 
the second-difference constraint at the point source by setting the rows of Q corresponding to 
the point-source location to zero. Figure 1./ shows a restoration of the same test image when 
the second difference constraint is not applied at the point source. A significant improvement 

is apparent. 


A direct extension of the method to 
two dimensional images is to minimize 
the Laplacian at each point. The Lapla- 
cian operator has a value at each pixel 
equal to four times the pixel value minus 
the values of the four immediate neigh- 
boring pixels. We use the subscript, L, 
to indicate the presence of the Laplacian 
constraint. As before, we set rows of the 
matrix Q to zero when the Laplacian con- 
straint is not appropriate (i.e. edges or 
point sources). 

The constraint need not be binary: 
we can vary the amount of constraint be- 
tween no constraint to full constraint for 
anv pixel, simply by multiplying the ap- 
propriate row in 0 by a constant factor 
running from 0 to 1. 

Another useful constraint is to mini- 
mize the difference of x from a trial solu- 
tion (i.e. minimize ||p - x|j). The solu- 
tion using Lagrangian multipliers is given 
by (Twomey 1963): 

X - (H t H + lt I)- l {H T b 4- 7tP) (6) 

where p is the trial solution, I is the iden- 
tity matrix, and 7 t is the reciprocal La- 
grangian multiplier. The subscript, t, will 
be used to identify the constraint as min- 
imization of the solution from a trial so- 
lution. Some possible choices for the trial 
solution, p, are a constant value (t.e. all 
zeros) or the blurred image itself. In ei- 
ther case, the ill-conditioned nature of H 
can be avoided and reasonable solutions 
obtained. 
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Multiple image constraints can be applied simultaneously: 

x = (H T H + 7 a Q T Q + • • • + 7 t IT l {H T b + 7 e P ) 

where a different value of 7 can be selected for each constraint. 

Selection of the reciprocal Lagragian multipliers can be done by visual inspection of the 
results for various values or by examination of the difference of blurred image and the solution 
re-convolved with the point spread function. This difference should have the same properties 
as the noise. 

2.3 Missing or Bad Data Values 

A problem occurs when trying to restore images with missing or bad data values ( i.e. cosmic 
ray hits or bad CCD columns). If these defects are not taken into account in the restoration, 
their bad values will propagate to a larger portion of the output solution. (To some extent, 
every point in the solution depends on all values in the blurred image.) 

One method of handling bad pixels is to attempt to repair them before restoration by 
interpolating from neighboring values. This approach is successful onl\ if the repair is accurate. 
An alternative method is to make no attempt at prior repair but handle them in the restoration 
process. In this approach, the restored image will have more data values than the blurred 
image, and the linear system is underdeterrmned and. therefore, singular (t.c. no direct inverse 
exists). To ignore defective pixels, we set the corresponding rows in matrix H to zero. 

This method of implementation (as opposed to removing row H creating a non-square un- 
derdetermined system) allows us to keep the matrix H square and decrease the complexity of 
implementation. Keeping H square in no way alleviates the problem of singularity. However, 
the method of constrained least squares solution does alleviate the problem of singularity and 
obtains reasonable solutions. 

3. RESULTS 

3.1 Ground Based image of QSO 2130+099 

Although the first example is of a ground base image, it illustrates a case that may come up 
frequently with Hubble Space Telescope (HST) data. We have a bright point source on a lower 
level diffuse source. The wings of the bright source makes it difficult to study the underlying 
diffuse structure. One question, which may prove crucial to the understanding of the origin of 
QSO’s, is: What kind of galaxy plays host to a QSO? The difficulty in answering this question 
is that the host galaxy appears as a faint fuzz around the bright QSO. 

Our approach to deconvolve the QSO image is to make the assumption that the center 
of the galaxy contains a point source. As described in section 2.2, we use the method of 
constrained least-squares, apply two constraints simultaneously, one involving smoothness in 
the restored image (eq. 5), the other involving the deviation from a trial solution (eq. 6 ). The 
smoothness constraint is appropriate for the host galaxy and background sky. It is empatically 
not appropriate for the nucleus (QSO), since that by definition is a point-source. This is where 
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Figure 2. Contour plots of Quasar 2130+099. Left: the original image. The 
maximum count (at the nucleus of Q2130+099) Is 21900 counts/pixel. Right: the 
restored image. The count level at the nucleus (the quasar) is now 324,000 
counts/pixel. Both images have a plate scale of 0.6 arcsec per pixel; thus 
Doth span a 1* x l 1 field. The two contour levels are at 280 and 320 counts 
per pixel; the average sky level is 254 counts/pixel. 



Figure 3. Cross-sectional plot of Q2130+099. The solid line shows the 
restored image. (The maximum count actually goes off-scale to 324,000 counts.) 
The pluses show the original (blurred image), while the dashed curve shows the 
result of convolving the restored image with the point-spread function. 





the algebraic approach is so useful: it allows us to have local control of the constraints We 
apply the smoothness constraint (minimize the Laplacian) and trial solution constraint (trial 
solution = sky background) to all the pixels in the image except the QSO. 

Our example is a CCD image of QSO 2130+099 obtained by Tim Heckman at the 4-meter 
telescope at CTIO. QSO 2130+099, also known as II Zwicky 136, is a relatively nearby QSO 
with a redshift of only 0.06. Figure 2-left shows a contour plot ot the observed image at its 
lowest count levels. Not only can we see the galaxy extending 0.5 arcmin across the sky, but we 
can see two protrusions from the nucleus that look like spiral "stumps’ if not full spiral arms. 
What we seek from deconvolution of Q2130+099 is not so much to enhance the resolution as 
it is to remove the veiling of the host galaxy by the QSO. Ideally, we would like to suck up 
all the flux from the QSO (nudues) into a single pixel, so that we can look at what is aroun 
it. Figure 2-right shows the contour plot of the restored image at the same contour levels as 
before. Now the spiral arms arms are more prominent and fully developed. 

Figure 3 shows a cross-sectional plot of the restored image. It. shows Q2130 099 for v, a 
it is- an exceedingly bright nucleus (324,000 counts) embedded in a galxy whose surface bright- 
ness falls off exponentially with increasing distance from the center, a brightness distribution 

typical of spirals. 


3.2 Wide Field Camera Image of R136 

Figure 4-left shows a Wide Field Camera image of R136 in the Large Magellanic Cloud. It 
shows a crowded field of stars embedded in an underlying halo resulting from the wings of t e 
HST point spread function. We could use the approach of the previous example: measuring 
the locations of all of the stars and applying a smoothness constraint at all locations except t e 
star locations. This, however, would be difficult for very close stars and very dim stars which 
are difficult to see in the image. Errors in the stars’ locations will result in a solution with 
artifacts. Another approach is to decrease the constraint at pixels with a higher probability 
of containing a star. A simple measure of the probability is the flux in the pixel. The larger 
the flux, the more probable that the pixel contains a star. In the solution, (figure 4-nght ) «e 
have used a weighted constraint that minimizes the norm (sum of the squares) of the difference 
of the solution from zero. The weight at each pixel was selected as the log of its value in the 
previous iteration divided by its value in the previous iteration. As the solution converges the 
constraint converges to the minimization of: 

Y^Xilog{xi) 

i 

Note that the solution is no longer Unear. Brighter stars are much sharper than dimmer stars. 
This can cause problems when performing photometry in the restored image. The accuracy 
of the photometry in the restored image is of major importance and will require additional 
investigation before we can make a reasonable assessment of the success of our restoration. 




Figure 4: Restoration of Wide Field Camera image of R136 

3.3 Wide Field Camera Images of Satvrn 

Figure 5 shows the restoration results for Wide Field Camera images of Saturn taken with 
three different wavelength filters. The images on the left are the unrestored, bias-subtracted 
and flat-fielded images. The images on the right show the results of the block iterative restora- 
tion algorithm using the minimization of the Laplacian (7^=0.001)and the minimization of the 
difference of each iteration from the previous iteration (7 ( =0.01) to constrain the solution. No 
attempt was made to repair bad pixels (e.g. cosmic rays) in the raw data. Instead, their loca- 
tions were manually flagged and treated as missing data as described in section 2.3. 

3.4 Goddard High Resolution Spectograph 

The last example is a one-dimensional spectrum taken by the Goddard High Resolution 
Spectrograph (GHRS). The GHRS has two square science apertures: a 2.0 x 2.0 arcsecond 
Large Science Aperture (LSA) and a 0.25 x 0.25 arcsecond Small Science Aperture (SSA). The 
apertures have the feature that they cut off the wings of the HST point spread function at the 
edges of the apertures. The result is that the HST spherical aberration causes almost no loss of 
resolution when the SSA is used but does cause a significant loss of light. Only approximately 
15 percent of the light from a point source centered in the SSA passes through the aperture. 
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Figure 5: Restorations (right) of bias subtracted and fiat-fielded WFC images of Saturn. Top 
- filter F439W. Middle - filter F547M, Bottom - filter F718M 
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Figure 8: Comparison of GHRS restored and unrestored LSA spectrum with a SSA spectrum 


Approximately 60 percent of the light passes through the LSA at the cost of a significant amount 
of resolution. 

A user of the GHRS must use the SSA (with a significant light loss) to obtain the besl 
resolving power unless deconvolution of LSA spectra can recover the resolution of SSA spectra. 
Figure 6 shows a spectrum of £ Persei observed through the LSA. Figure 7 shows the results 
of the block iterative restoration (solid line) compared to a SSA spectrum with a much smaller 
signal to noise (dots). Not only does the restoration separate the P II and 0 I lines at approxi- 
mately 1304.8 angstroms (Also see figure 8) but it correctly restores the profile of the Si II line 

at 1304.3 angstroms. 
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